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The relaxation dynamics of the one-dimensional totally asymmetric simple exclusion process on 
a ring is considered in the case of step initial condition. Analyzing the time evolution of the 
local particle densities and currents by the Bet he ansatz method, we examine their full relaxation 
dynamics. As a result, we observe peculiar behaviors, such as the emergence of a ripple in the 
density profile and the existence of the excessive particle currents. Moreover, by making a finite-size 
scaling analysis of the asymptotic amplitudes of the local densities and currents, we find the scaling 
exponents with respect to the total number of sites to be —3/2 and —1 respectively. 

PACS numbers: 05.60.Cd, 02.50.Ey, 75.10.Pq 



Introduction.— Nonequilibrium statistical mechanics 
has been a hot topic for the last twenty years. Although 
a large amount of research has been done, the theory 
beyond the linear response regime is still far from being 
complete. One strategy to gain insights into nonequilib- 
rium statistical mechanics is to examine special models 
in which exact methods can be employed, and to extract 
as much information as possible from them. 

Among them, the asymmetric simple exclusion pro- 
cess (ASEP) [IH| is one of the well-studied paradigms 
in nonequilibrium statistical mechanics. The ASEP is a 
continuous time Markov process describing the diffusion 
of classical particles on a one-dimensional lattice. Each 
particle is viewed as a biased random walker that obeys 
an exclusion principle: each lattice site can not be occu- 
pied by more than one particle. This simple model ap- 
pears in many different contexts. It originally appeared 
to describe the dynamics of ribosomes along RNA 
and have modern applications to pedestrian and traffic 
flows 0, , transport of quantum dots [8] . It can also be 
regarded as a discrete analogue of the KPZ equation 0|, 
which describes the surface growth phenomena, observed 
in recent experiments on electroconvection fioj . 

The ASEP is rather a toy, but plenty of exact results 
have been obtained and still continuously provides valu- 
able insights into nonequilibrium statistical mechanics. 
For instance, as for the steady state properties, the ma- 
trix product ansatz was used to reveal various interesting 
phenomena such as the boundary induced phase transi- 
tions pj]-[m. On the other hand, as for the dynami- 
cal properties, the asymptotic behavior of the relaxation 
process was exactly investigated by the Bethe ansatz 
[2l|,[28j, and shown that it is governed by the KPZ univer- 
sality class. In addition, the current fluctuation, which 
is known to be described by the largest eigenvalue dis- 
tribution of random matrices, is recently studied by use 
of random matrix theory and the Bethe ansatz |22|-[28|. 
Despite these remarkable developments, it is still a chal- 
lenging problem to examine the exact relaxation dynam- 
ics beyond the first principles Monte Carlo simulations 
which requires many samples to obtain reliable results. 



In this article, we study the exact dynamics of the to- 
tally asymmetric simple exclusion process (TASEP) on a 
periodic ring, which is a special case of the ASEP where 
particles are allowed to move only in one direction. Start- 
ing from the step initial condition where the half of the 
system is consecutively occupied by the particles and the 
other half is empty, we investigate the relaxation process 
to the steady state by use of the Bethe ansatz method. 
By quantitively evaluating local particle densities and 
currents, we observe interesting behaviors peculiar to the 
step initial condition such as the emergence of a ripple in 
the density profile and excessive flow of currents. This is 
the first study to evaluate the full relaxation dynamics of 
the ASEP by using the Bethe ansatz method. Moreover, 
by making a finite-size scaling analysis of the asymptotic 
amplitudes of the local densities and currents, we find 
the scaling exponents with respect to the total number 
of sites to be —3/2 and —1 respectively. The results 
are consequences of the advantages of the Bethe ansatz 
method which allows us to conduct the finite-size scaling 
of the amplitudes completely separate from the exponen- 
tial parts determining the relaxation times. 

Definition of the TASEP.- We consider the TASEP 
on a periodic ring with M sites and N particles. Each 
site can be occupied by at most one particle (see Fig. [I]). 
The dynamical rules of the TASEP is as follows: during 
the time interval dt, a particle at a site j jumps to j + 1th 
site with probability dt, if the j + 1th site is empty. A 
Boolean variable T{ is associated to every site j to indicate 
whether a particle is present (t; = 1) or not (ji = 0). The 
probability of being in the (normalized) state |n, . . . , tm) 
is denoted as Pt(ri, . . . , tm). The time evolution of the 
state vector \ip(t)) = J2 Ti =o,i p t( T i: • • • , r M )\n : . . . , r M ) 
is subject to the master equation 
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FIG. 1. The step initial condition of the TASEP. TASEP is a 
special case of the ASEP where particles are allowed to move 
only in one direction. 
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It can be shown that the above states are eigenstates of 
the Markov matrix (|2]), if the spectral parameters {u} = 
{^1,^2, ...,un} satisfy the Bethe ansatz equation 



N 



(i 



2 ) 



-M -2N 



(-if- 1 !!^ 



(6) 



Here the Markov matrix Ai of the TASEP is defined by 
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where := (<jj ± icr y -)/2 and a x -' V)Z is the Pauli ma- 
trices acting on the jth site. (Note here that we inter- 
pret that cr||ri,...,r M ) = (-l) Tj \n, . . . , r M )). The dy- 
namical rules of the TASEP are encoded in this Markov 
matrix. To describe the relaxation dynamics, we must 
evaluate the eigenvalues as well as the eigenstates of the 
Markov matrix Ai. The algebraic Bethe ansatz is one of 
the most useful methods to achieve this. 

Algebraic Bethe Ansatz for Dynamics.— The Markov 
matrix A4 ([2]) can be interpreted as a (non-Hermitian) 
Hamiltonian of a one-dimensional quantum spin system 
which is known to be integrable. Therefore various exact 
methods developed in the study of integrable systems are 
also available to analyze the dynamics of the TASEP. 

To describe the relaxation dynamics by the Master 
equation (pQ) , one must calculate both the eigenvalues and 
eigenstates of (pQ) because the time evolution of the phys- 
ical quantities are, in general, given by the form factors of 
quantum operators (see ©). To this end, we employ the 
algebraic Bethe ansatz 0, [30| which makes it possible 
to calculate the correlation functions and form factors of 
quantum integrable systems. 

What plays a fundamental role is the L-operator acting 
on the jth site 

L(j\u) = ussj + n(ul — u~ 1 Sj) + a~ a+ + cr+aj , (3) 

where Sj = (1 + cr|)/2, nj = (1 — cr|)/2 is the projection 
operator onto the empty and filled states at jth site, re- 
spectively. The monodromy matrix is defined by a prod- 
uct of L-operators: 
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The arbitrary TV-particle state (resp. its dual (^|) 
(not normalized) is constructed by a multiple action of B 
(resp. C) operator on the vacuum state \Q) := |0, . . . , 0) 



for k = 1, 2, • • • , N. Then the eigenvalues of the Markov 
matrix are given by 
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Now we formulate the relaxation dynamics of the 
TASEP within the algebraic Bethe ansatz. The time evo- 
lution of the expectation value for the physical quantity 
A starting from an initial state \In) is defined as 



(A) t = (S N \Ae Mt \I N ) 



(8) 



Here \Sn) (resp. (Sn\) is the A^-particle steady state 
(resp. its dual) which is simply constructed by the su- 
perposition of all configurations with equal probabilities. 
We shall study the relaxation dynamics starting from 
the step initial condition (see Fig. [1]) where the half of 
the system is consecutively occupied by the particles and 
the other half is empty. The (normalized) initial state 
\I N ) = | 1 I _^,0 1 _^0) is given by \I N ) = B(1) N \Q) 

N M-N 

(see [31] for the XXZ spin chain). (Note that this ini- 
tial state is not the eigenstate of the Markov matrix). 
Thus the local densities = (1 — s*)t and currents 

(ji)t = ((1 - Si)s i+1 )t are respectively given by 
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Here the sum in the above is performed over all states 
except for the steady state (denoted by the index a). 

The remaining problem is the evaluation of the norms 
and the scalar products of the Bethe vector, and the form 
factors of the local operators. These quantities can be 
calculated in the framework of the algebraic Bethe ansatz 
(cf. [30j] ) . For the norm of the Bethe vector one 
obtains 
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with N x N matrix 

N — 1 + (M - N + l)u~ 2 

Qjl = ~2 — ^ " (1 " Sjl). (11) 
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On the other hand, the overlap between the initial state 
and the arbitrary Bethe vector is reduced to the following 
simple form: 
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Finally the form factor for the local operators is explicitly 
given by 



N 



N 



{S N \ Sl • • • Si+fc _!|V) = ~ Uj 2 )^- 1 ]^[ Uf +1 

i=i j=i 

N>l>n>l U l Un 

where the N x N matrix V is written as 
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Note that the overlap between the steady state and the 
Bethe vector (Sn\iJj) is obtained by setting i = = 
in (p~3|) . Inserting all of the above formula into (|9]), one 
can calculate the time evolution of the expectation values 
of the local particle densities and currents. 

Results.— We use the algorithm in [l8| to compute the 
roots of the Bethe ansatz equation (|6)) for the arbitrary 
states as many as possible, and insert into the formu- 
las (|9j), ([TO]) . (fT2|) and §F5§ to evaluate the local densi- 
ties and currents. We shall give a remark here for the 
algorithm in [18], which conjectures that one choice of 
a monotonous function gives one Bethe roots. We find 
this one-to-one correspondence does not hold anymore 
for M > 10 at half-filling. There may be two or more 
or zero Bethe roots for one choice of quantum numbers. 
However, these irregular Bethe roots give highly excited 
states. The missing roots are about 1% and does not 
not seem to contribute much to the physical quantities, 
as can be seen from the fact that the density profile of 
the initial state is almost completely realized. Here we 
demonstrate the case for 20 sites and 10 particles (half- 
filling). 

The top panels of Figs. [2]and[3]show the time evolution 
of the expectation values of the local particle densities. In 
the beginning, the behavior of the local densities matches 
with our intuition: as the site is closer to the 10th site 
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FIG. 2. (Color) The time evolution of local particle densities 
(rn)t and currents (ji)t for 20 sites, 10 particles. 



where the top particle of the line occupies, the local den- 
sity of a site converges faster to 1 /2 which is the density 
of the steady state. This is because the particles begin to 
flow around the 10th site. However, some time later, the 
density on the 1st site becomes to decay faster than the 
4th site, for example. This behavior can be intuitively 
explained: No particle comes into the 1st site until the 
top particle which initially occupies the 10th site reaches 
the 20th site. More precisely, the decay rate of the ex- 
pectation value of the particle densities is described by 
the continuity equation: 
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For k is small (say k = 1), the current does not flow 
for a while. After some time, (ji) becomes visible, but 
O20 )(= O'o)) is still very small: (j 20 ) <C (ji) (see the 
bottom panel of the Fig. [2]). Thus, the density of the 1st 
site (ni) decreases in a rapid way after when the current 
(ji) begins to flow. In contrast, for the sites where the 
middle particles initially occupy (say k = 4), the current 
0*3 ) begins to flow shortly after (74) has started to flow. 
The existence of the entering flow (js) moderates the 
decrease of the density of the 4th site (724). After some 
time, its value becomes larger than that of the 1st site 
where no particle enters until the top particle reaches. 
In this way, we observe the emergence of a ripple in the 
density profile. 

The bottom panels of Figs. [2] and [3] are the relaxation 
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FIG. 3. (Color) The time evolution of density profile (rii)t 
and current profile {ji)t for 20 sites, 10 particles. 



process of the local currents. First, note that the current 
profile is symmetric with respect to the 10th site. This 
reflects the fact that we are considering the half-filling 
case, and a system with particles moving from left to 
right with initially occupying 1, 2, • • • , 10th sites can be 
regarded as a system with holes moving from right to left 
with initially occupying 11, 12, • • • , 20th sites. 

The current of the 10th site is 1 at t = since the 
11th site is vacant, and monotonically decreases to the 
steady state current. The currents around the 10th site 
show a rapid increase from to reach peak, and then 
decrease. These mean that one observes the excessive 
flow of currents around the 10th site which the top par- 
ticle occupies in the initial state. As the site is more 
distant from the 10th site, the increase slope of the cur- 
rent gets flatter, and finally the peak disappears, i.e., 
the current no longer flows excessively and just shows 
monotonic increase. The behavior that the currents have 
peak in some regime between the 1st and 10th sites can 
be regarded as the combined effect of the monotonic de- 
crease represented by the 10th site where the top particle 
immediately starts to move and the monotonic increase 
represented by the 1st site where the last particle has to 
wait for a while to start moving. 

Finite-size scaling.— We conduct the finite-size scaling 
analysis at half-filling. For large times, local densities 
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FIG. 4. (Color online) The lowest excited energy .Mist plotted 
against the total number of sites M for M — 2 k , k — 1, • • • ,10. 



and currents behave as 
N 
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where A4i s t is the non-zero eigenvalue of the Markov ma- 
trix M with largest real part. First, as a check, we 
perform the scaling analysis of A^ist- Fig. |4] shows 
ln(— A^i s t) vs InM. The simplest fitting from M = 
256, 512, 1024 gives \n(-Mi st ) = 1.89793 - 1.503511nM, 
which implies the KPZ scaling (see pjl [20| for the thor- 
ough analysis of the open as well as the periodic boundary 
conditions of ASEP). 

Next, we make the scaling analysis of the asymp- 
totic amplitudes A(rik) and A(j k ) of the local densi- 
ties (rik)t an d currents (jk)t- We take the amplitudes 
at the N-th site A(un) and A(jn) as a representa- 
tive. The top and bottom panel of Fig. [5] shows 
ln(A(n/v)) vs InM and \n{A{jisr)) vs InM, respectively. 
The fitting from M = 256,512,1024 gives \n(A(n N )) = 
1.937981 - 1.49331871nM and ln(A(j N )) = 0.854409 - 
0.99685791nM. Confining to larger sites for fitting, the 
coefficient associated with InM approaches to —3/2 and 
-1, which indicates A(un) oc M -3 / 2 and A(j N ) oc M _1 
respectively. 

Conclusions.— In this article, we have studied the ex- 
act relaxation dynamics of the TASEP. We examined the 
behavior of the local densities and currents by use of the 
algebraic Bethe ansatz method. By quantitive evalua- 
tion, we see the convergence of the density and current 
profiles to the steady state densities and currents. Fur- 
thermore, we observe interesting behaviors such as the 
emergence of a ripple in the density profile and exces- 
sive flow of currents peculiar to the step initial condition. 
Moreover, by making a finite-size scaling analysis, we de- 
termine the scaling exponents of the asymptotic ampli- 
tudes of the local densities and currents with respect to 
the total number of sites, which are found to be —3/2 
and —1 respectively. 



5 



2 

-8 
-10 

































































• 


m 









































































































































































































































































































1 2 3 4 5 6 7 







Li' 


VI 




V 















































































































































































































































1 2 3 4 5 6 7 

InM 



The authors would like to thank C. Arita for useful 
discussions. The present research is partially supported 
by Grant-in-Aid for Young Scientists (B) No. 21740285. 
J.S. is supported by JSPS. 



FIG. 5. (Color online) The asymptotic amplitudes A(un) and 
A(jn) of the local density (rijv)t and current (jN)t plotted 
against the total number of sites M for M = 2 k , k = 1, • • • ,10. 
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